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Abstract 

A new non- Fermi liquid state of quasi-one-dimensional conductors is suggested in which electronic 
system exists in a form of collection of bounded Luttinger liquids stabilized by impurities. This 
state is shown to be stable towards inter-chain electron hopping at low temperatures. Electronic 
spectrum of the system contains zero modes and collective excitations of the bounded Luttinger 
liquids in the segments between impurities. Zero modes give rise to randomly distributed localized 
electronic levels, and long-range interaction generates the Coulomb gap in the density of states at 
the Fermi energy. Mechanism of conductivity at low temperatures is phonon-assisted hopping via 
zero-mode states. At higher voltages the excitations of Luttinger liquid are involved into electron 
transport, and conductivity obeys power-law dependence on voltage. The results are compared 
with recent experimental data for NbSe3 and TaS3 crystals. 

PACS numbers: 71.10.Pm, 71.10.Hf, 71.27.+a, 71.45.Lr 



I. INTRODUCTION 



The concept of single-electron quasiparticles is one of the central ideas of solid state 
physics. Indeed, systems of interacting electrons in metals and semiconductors are well de- 
scribed by Landau's Fermi-liquid picture in which interaction modifies free electrons making 
them quasiparticles that in many respects behave like non-interacting electrons. Basic elec- 
tronic properties of many solids, including phase transitions into symmetry-breaking states 
like superconducting or charge- density wave (CDW) state, are well described within this 
concept. Therefore, possibility to observe non- Fermi liquid behavior in conducting materi- 
als attracts much interest. The concept of Luttinger liquid is an alternative to the Fermi 
liquid elaborated for one-dimensional (ID) electronic systems. It was found that in ID 
electronic systems the Fermi-liquid picture breaks down even in case of arbitrarily weak 
interaction. Single-electron quasiparticles cannot exist in ID metals, and electrons form the 
Luttinger liquid in which the only low energy excitations turn out to be charge and spin 
collective modes with the sound-like spectrum. Dynamical independence of these modes 
gives rise to a spin-charge separation in ID systems. Furthermore, correlation functions at 
large distances and times decay as a power law with interaction dependent exponents (for a 
review see Ref. 1,2,3). 

One may expect to observe formation of Luttinger liquid in quasi-lD conductors, i. e., 
highly anisotropic 3D conductors with chain-like structure. However, according to present 
theoretical point of view, the formation of Luttinger liquid in quasi-lD conductors at low 
enough temperatures is problematic because of the instability towards 3D behavior in the 
presence of arbitrarily small inter-chain hopping 4 ' 5,6 ' 7,8 ' 9 . So inter-chain hopping induces 
a crossover to 3D behavior at low energies, while Luttinger liquid behavior can survive 
only at high enough energy scale where it is not affected by 3D coupling. In contrast to 
inter-chain hopping, the Coulomb interaction between electrons on different chains does not 
destroy the Luttinger liquid state, it merely modifies the electronic spectrum of the Luttinger 
liquid 10,11,12,13 . 

Numerous experimental studies of both organic and inorganic quasi-lD conductors at low 
temperatures demonstrate, commonly, Fermi-liquid metallic behavior and/or transitions to 
broken-symmetry states, which are described in terms of Fermi-liquid ideas (for a review 
see Refs. 14,15). For instance, the most studied inorganic quasi-lD metals (e. g., blue 
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bronze K0.3M0O3, TaS 3 , (TaSe 4 ) 2 I etc.) undergo the Peierls transition from metallic to 
a semiconducting CDW state, or to semimetallic CDW state like NbSe 3 . Typically, this 
transitions occur in the temperature range 50 -j- 250 K. Some quasi- ID materials, e. g., 
TaSe 3 do not undergo the Peierls transition and remain in the normal metallic state at low 
temperatures as well. 

However, a transition from room-temperature metallic behavior to nonmetallic one ac- 
companied by disappearance of the CDW state at temperatures below 50 100 K was 
detected in recent experimental studies of temperature and field dependence of conductivity 
of TaS 3 and NbSe 3 nanoscale-sized crystals 16 ' 17 ' 18 . The low temperature non-metallic state 
was characterized by power law dependencies of conductivity on voltage and temperature like 
that expected in Luttinger liquid, or by more strong temperature dependence corresponding 
to the variable-range hopping 18 . Resembling dependencies of conductivity were reported 
also in focused-ion beam processed or doped relatively thick NbSe 3 crystals 19 . Transition 
to non-metallic conductivity at low temperatures was observed earlier in fragmented into 
small pieces and then sintered NbSe 3 crystals 14 . Hopping conductivity was also found out in 
heavily doped by iron bulk NbSe 3 crystals 20 in dirty quasi-lD conductors KCP and organic 
TCNQ-based metals 21 , while pure materials are known to undergo the Peierls transition to 
the CDW state. 

In order to account for such behavior a possibility of stabilization of the Luttinger state by 
defects in quasi-lD metals was put forward recently in the letter of one of the authors 22 . Such 
a possibility is based upon finding made in Refs. 23,24,25 that a local impurity potential in 
Luttinger liquid acts, at low energies, as effectively infinite barrier. This leads to a discrete 
spectrum of collective charge and spin density fluctuations, so that the electronic system 
can be considered as a collection of bounded Luttinger liquids. In this case, at temperatures 
below the minimum excitation energy of the collective modes, weak inter-chain hopping can 
be considered as small perturbation, and, therefore, hopping does not destroy the Luttinger 
liquid state. 

However, in order to make comparison with the experimental data and to find further 
evidences that the low-temperature non-metallic state observed in quasi-lD conductors can, 
actually, be interpreted as a collection of impurity-induced bounded Luttinger liquids, cal- 
culation of conductivity in the latter state is needed. Such calculations are the main goal 
of this paper. We calculate conductivity in a collection of weakly coupled metallic chains 
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with defects in the limit of low temperatures, considering both the case of the contact in- 
teraction used in the standard Tomonaga-Luttinger model 1 ' 3 and the case of more physical 
long-ranged Coulomb potential. 

The paper is organized as follows. In Sec. II we present main equations and formulate the 
problem in terms of bosonization technique generalized to the case of long-range Coulomb 
interaction between electrons in multi-chain problem in presence of impurities. In Sec. 
Ill we study electronic structure of quasi-lD conductors. We consider, first, interacting 
electrons neglecting inter-chain hopping and find the solution for the phase operators in 
the limit of strong impurity potential ignoring, at the first step, the Coulomb interaction of 
electrons at different chains. Modifications induced by the long-range nature of the Coulomb 
interaction are discussed in the next subsection. In the last subsection we show that at low 
temperatures the inter-chain hopping can be considered as small perturbation, provided the 
hopping matrix element is small enough. Sec. IV is devoted to transport properties of 
the system. We show that in the limit of low voltages conductivity is described in terms 
of variable-range hopping, while at large voltages it turns into conductivity described by 
a power-law dependence on voltage. And, finally, in Sec. V we discuss relation of the 
theoretical results to experimental data, and make conclusions. In this paper we use units 
with Planck and Boltzmann constants equal to unity, H — 1, ks — 1. 



II. GENERAL EQUATIONS 

First of all we start with the Tomonaga-Luttinger model ignoring inter-chain hopping 
integral, t±, in order to formulate the problem in the zero approximation in t±. Electronic 
operators for right (r = +1) and left (r = —1) moving electrons with spin s on chain n are 
given in terms of phase fields as (see Ref.1,3,11) 

^„, s (r, x) = lim - 7 = Vn , ryS e- iA ^ x ' n \ (1) 
A n ,r = -j=[9 p (x, n) - r<S> p (x, n) + s(O a (x, n) - r$ a (x, n))]. 

Phase fields Q u (x,n) are related to charge {y = p) and spin {y = a) densities, while fields 
Q v (x) are related to the momentum operators IIj, = (l/n)d x Q u canonically conjugate to 
Further, rj r jS are Majorana ("real") Fermionic operators that assure proper anticommutation 
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relations between electronic operators with different spin s and chirality r, and the cut off 
length a ~ 1//cf is assumed to be of the order of interatomic distance. 

The Hamiltonian of the system includes three terms, the Hamiltonian of the interacting 
electrons, the impurity scattering term, and the inter-chain hopping Hamiltonian. The first 
term containing kinetic and potential energies of the electrons is described by the Tomonaga- 
Luttinger Hamiltonian generalized for the multi-chain case with long-range interaction that 
couples electrons on different chains 1,3 ' 11 . This Hamiltonian does ignores the backscattering 
terms in the electronic coupling, i. e., it neglects components of interaction potential with 
q\\ ~ 2kp. In the bozonized form the Hamiltonian reads 

where d 2 is the area per single chain, and integration over is taken within the first 
Brillouin zone 

**(cu,?||) = (l + — V^, (3) 
K v is the standard Luttinger liquid parameter describing the strength of the interaction. 
K v determines velocities of the charge (is = p) and spin (y = a) modes, v v = vf/K v . We 
study the case of spin-independent interaction, therefore, in our case, the coupling constant 
for spin channel is equal to zero, g a = 0, and K a = 1. For charge channel <? p (q) coincides 
with the matrix element measuring the strength of the forward scattering due to interaction 
between electrons on the same or on different branches of the electronic spectrum. The exact 
form of <7 p (q) depends on crystalline and electronic structure of a crystal, but in the long- 
wavelength limit, gdc 1, one can describe interaction by the standard Fourier transform of 
the Coulomb interaction potential. 

. 2 2 

2 2 ^ vl 111 / a \ 

fr = dV ± + g f) ' vp = VF + &^{ up = q " vp) () 

where lo v \ is the plasma frequency and uo p is the frequency of the charge mode. 
In the coordinate representation 

2e 2 

K; 2 (x, n) = 5 n0 5(x) + (5) 

iivfvx + n z d A 

where the last term is induced by the Coulomb interaction. Its contribution to the Hamil- 
tonian (2) in the coordinate representation reads 

EE/ dxdx ' | n) f — d x '®„(x', n') 1 (6) 

^=P^n,n' J [ 7T 2 ^(X - X>) 2 + (ll- Il') 2 d 2 J 



Since the operator of the particle density is given by expression p = — (y/2/n)d x $(x), this 
term has rather transparent physical meaning. 

The second part of the total Hamiltonian is 2k F impurity backscattering term 1 ' 2 ' 3 

H imp = -Y, I dx ~ cos (V2$ p + 2k F x) cos (y/2Q ff (x)), (7) 

n J i 7101 

where V(x) is potential of impurity centered in position i. In calculations below we assume 
V(x) ~ V X5(x), where V and A ~ a are amplitude and radius of the scattering potential, 
respectively. 

When inter-chain electron hopping is ignored, equation of motion for Heisenberg phase 
operators <3> p and 6 P can be obtained in standard way from the Hamiltonian H + H imp 
given by Eqs. (2) and (7): 

d t e p (t,x,n) =v F Y.J dx'K- 2 (x - x',n - n')d x/ ^\t,x',n')t, (8) 

n' 

<9 2 $ p (t, x,n)-v F J2j dx'd x K~\x -x',n- ri)d x ,$ p (t, x', n') 

n' 

= ^ VodvF Y,6(x- Xi ) sin (y/2$ p + 2k F x) cos ( V2$ a ) . (9) 

a i 

Kane and Fischer 23 found that the backscattering impurity potential for repulsive poten- 
tial (K p < 1) flows to infinity under scaling. This conclusion was made for the case of single 
chain and a contact interaction potential. 

The results of Refs. 23,24,25 were generalized by Fabrizio and Gogolin 26 to the case of 
many impurities. It was shown that the impurity potential can be considered as effectively 
infinite provided that the mean distance, I, between impurities satisfies the condition 

1 /D\ 2/{1 ~ Kp) 

J »r F k) ■ < 10 > 

where D is the bandwidth. We assume that the impurity potential is of atomic scale, Vq ^ D, 
and the interaction between electrons is not too weak, (i. e., K p is not too close to 1). Then 
condition (10) is satisfied for I ^> l/k F , so the limit of strong impurity potential should be 
a good approximation in a wide range of impurity concentrations. 

Similar conclusions can be derived for the multiple-chain system when interaction between 
electrons on different chains is taken into account. Indeed, calculations in Ref. 23 are based 
upon dominant contribution to Euclidian action expressed in terms of the value of the field 
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$ p at the impurity site by means of solution of equation of motion in imaginary time. One 
can act, similarly in case of multiple-chain system. We express the solution of imaginary- 
time version of Eq. (9) for Fourier component of $ p in terms of its value at the impurity 
site, $0; an d find 

$ ( f d 2 dq±dq\\ 1 

[J (2vr) 3 u£ + g}i 

In the spirit of renormalization group approach the last factor in Eq. (11) must be calculated 
in the limit of small energies. Since at small u> this integral diverges at small q\\, we use 
v 2 = Vp/K 2 at 5|| =0. Then we obtain 

d 2 dq±dq\\ 1 K 



/ 



(2vr)3 ou 2 n + q 2 v 2 2v F u' 
where 

v lis ( n^\ I dq±K p ((i ± ,q ll = 0) 
K = (K p (q ± , qil = 0)) q± = fd^T ' 

Then using (11) one finds the action 

S = E H + 4 V > 2 (^, q)) = E 2 = v*&0). (12) 



This equation is similar to the action in Ref. 23 that was used to derive the renormalization 
equation for the interaction parameter, the only difference being that in case of long-range 
interaction action (12) contains the averaged value K of the interaction parameter instead 
of parameter K p itself. Thus impurity potential flows to infinity under scaling if K < 1. 

So at energies close to the Fermi energy the limit of strong impurity potential can be used. 
In the main approximation this limit corresponds to boundary conditions at the impurity 
site 

V2$ p + 2k F Xi = nil, V2$ a = mir, (13) 

where n and m are integers. Furthermore, n + m must be an even integer in order to ensure 
the minimum value of the impurity Hamiltonian (7). This connection violates independence 
of spin and charge modes similar to violation of spin-charge separation in bounded Luttinger 
liquid discussed by Eggert et al. 27 ' 28 . 

The last term to be added to the total Hamiltonian describes next-neighbor inter-chain 
hopping in the standard way 



H± = t± E / dx ^,s, m (x)A,sA x ) + h.c. 



E / dx %t±Vr ' 8 ' nVr ' a ' m [sin(A r , m - A,n) + sm(A, m - A. r , a + 2vr-k F x)\ , (14) 

n tj t* c 7TG! 



,s,m 

m,n,r,s ' 

where indices n and m denote the chain numbers related to the nearest neighbors. 



III. ELECTRONIC STRUCTURE AT LOW TEMPERATURES 

We consider first the case of interacting electrons neglecting inter-chain hopping and find 
the solution for the phase operators in the limit of strong impurity potential. According 
to discussion in the previous section this limit should be a good approximation in a wide 
range of impurity concentrations. In this approximation the system breaks up into a set 
of independent segments. At the beginning we ignore the long-range Coulomb interaction 
and consider the case of contact interaction as in the standard Tomonaga-Luttinger model, 
and then discuss modifications induced by long-range interaction. At the end of the section 
we will show that at low temperatures the inter-chain hopping does not produce qualitative 
modification of the electronic structure based on the concept of Luttinger liquid. 

We begin with solutions for phase operators in case of contact interaction, taking into 
account that 9^ and Q u must obey the commutation relations 1 ' 2 ' 3 ensuring anticommutation 
of electronic operators (1). Using then the analogy of H in Eq.(2) with the Hamiltonian 
of an elastic string strongly pinned at impurity sites, and taking into account commutation 
relations, we can write down solutions for the phase operators in the region between impurity 
positions at x = Xi and x i+i as 

*u = ^p- nAN ^ - ^ AN uj + £ J^{b nvi e-^ + b+ ut eS) sin (15) 

1 j<i n=l » 71 I" 1 



9„ = 0„ + ^ m> Nv + % E . -—(br^e-^ - b+ vi e^) cos — - , (16) 

where x = x — Xi, U = Xi+i — x-i. 

Solutions (15-16) consist of two parts. The terms with summation present a general 
solution with zero boundary conditions, they describe excitations. Excitation spectra of the 
eigenmodes are u nui = mrv l/ /l i . 

The first terms present the particular solution with boundary conditions (13), they de- 
scribe zero modes. In the zero modes 

„ _ «f _ v v (A% + AN lt ) (ANjj - AA it ) 
^ = A^ = AV AiV -- 71 ' m ~ 71 ' 



AiVjj (AJVjj) is the number of extra electrons with spin up (down) in the region between 
impurities number i and + V2S(f>i is equal to the modulo 2n residue of 2kpk, and, finally 
9 U i is the phase canonically conjugate to AN ui obeying commutation relations [9 U i, AN ui ] = i. 

Note that inside the segments between the impurities the expressions for the phase fields 
between the impurity sites turn out to be similar to those found for bounded ID Luttinger 
liquid with open boundary conditions at the sample ends 26 ' 27,28 . The main difference is 
that Majorana operators i] s in Eq.(l) are the same for electrons moving right and left, and 
that Eq.(15) contain the summation over j < i that insure proper commutation relations 
between the electron operators related to different segments. Thus we conclude that the 
system breaks up into a set of independent segments described as bounded Luttinger liquids 
with main properties similar to those discussed in Refs. 26,27,28. In particular, as long as 
eigenvalues of AN^ and AN^ are integers we find that AN pi = n pi /\^2 and AN ai = n ai /y/2 
are not independent, because n pi + n ai must be an even number. Such a limitation ensures 
the minimum value of the impurity Hamiltonian (7). This also means violation of spin- 
charge separation due to zero modes as discussed by Eggert et al. 27 ' 28 for bounded Luttinger 
liquids. 

Hamiltonian in the region between impurity positions at x = Xi and x i+ i can be presented 

as 

H = J2( e u + J2 ^nvbtnKn) + constant (17) 

v n 

where contribution of zero modes has a form 

= VNpjS^j - 7rAiV p ) 2 e = v F (7iAN a ) 2 
'2~l , '2~l ; 

One can see that energies of the zero-mode states depend on number of extra electrons 
at a segment confined by impurities. 

From the constitutive relation (1) one can see immediately that temporal dependence 
(16) of zero-mode part of 6 should result in local energy levels: 

(<J0i - nANpi)v Np + TiAN ai v Na ^ 



Vv,s(0 °c exp 



-^=(e niP + se ni(T ) 



oc exp 



-t 



V2h 

Strictly speaking, the single-electron energy levels can be found from poles of single-particle 
Green's functions. Its calculation is similar to that in Ref. 28. As long as we study the 
low-temperature behavior of the system, we need explicit expressions for Green's function 
at T <^ u ^. It is presented in Appendix. 
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A. Zero modes 



Now we concentrate on zero modes. Consider the zero-mode state (n p , n a ) with extra 
charge n p and extra spin n a s in one of the segments. Such a state contains (n p + n a )/2 extra 
electrons with spin s and (n p — n CT )/2 electrons with spin — s. Energy (18) for formation of 
the zero-mode state (n p ,n a ) can be conveniently presented in a form (in this subsection, for 
brevity, we drop the index related to the number of a segment). 

_ u 0p ^-n p ) 2 _ iv 0a nl _irv F 

where we introduced a random factor £ = a/250/7t, |£| < 1. Since factor £ is related to 2k F l 
and, hence, depends on positions of impurities confining the segment under consideration, 
it is different for different segments. 

The low-lying frequencies of the single-particle Green's function G + ~ and G~ + introduced 
by Keldysh (see Ref. 29) describe, respectively, energy levels of electron and hole states. So 
according to Eqs. (A2) and (A3) the energy levels induced by zero-mode states (n p ,n a ) are 

£n p n a = ^ [uQ P n p + uj 0a n a s - UJ 0p {£ ± f )] (20) 

where signs ± are related to electron and hole states, respectively, £o = \ [1 + (K p /K a ) 2 ]. 
(We remind that we consider spin-independent interaction, so in our case K a = 1). 

As long as typical values of £ and £ are of the order unity, a characteristic energy for 
formation of a zero-mode state (19) and characteristic energy scale for position of a zero- 
mode level (20) can be estimated as 

71 Vp 



u; 



with I = 1/N imp , where N imp is the impurity density per single chain. If impurities are 
distributed randomly, then probability to find a segment of length / is given by the Poisson 
distribution, w(l) = Ni mp exp (— Ni mp l) . Then, for a given segment, probability to find a 
value of frequency u that is much smaller, than ZJ , is exponentially small. Hence, at low 
temperatures, T vf/1, only energies of modes with (n p ,n a ) equal to (0,0), (1,±1), and 
to (—1, ±1) can be of the order of temperature. Other frequencies except mentioned above 
are at least by u larger. Furthermore, for a given value of £ not more then two states can 
have energies near the Fermi level simultaneously. 
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So in the statistically averaged zero-mode part of the Green's function presented in Ap- 
pendix we can keep only terms related to the states mentioned above: 

(e 7Ti ) e — = _( e *(9i,i^a-ei,i*a) + e i(gi,-ix„-ei,-it )) ni + 

_(gi(?-l,l*a- e -l,lta) + gi(9-l,-ia ! »- e -l,-lt»)) n _ 1 + g<(«),0X.- e 0,0t«)( 1 _ ^ _ ( 2 l) 

with q n na = + sn <r T 1 — 0- Distribution functions for electrons at zero-mode level 
are given by expressions 

H±1 = l + Iexp)^)]- (22) 
Note that the distribution functions must not coincide necessarily with Fermi-Dirac ones 
because electrons are correlated, and factor 1/2 in the denominator of Eq. (22) reflects 
presence of two spin states. 

Thus zero modes form the system of local electronic levels near the Fermi energy. The 
energies of these levels are determined by random factors |£| < 1, and density of these states 
at the Fermi energy is finite in case of short-range interaction. 

B. Modifications of zero modes by long-range interaction 

Zero modes in case of long-range Coulomb interaction can be found from equations of 
motion (8-9) with time derivative equal to zero, 

d x e?\x,n) = 0, (23) 
Q^(x,n) =0 V + Y, J dx'v F K~ 2 (x-x',n-n')d x ^(t,x',n% (24) 

n' 

with boundary conditions \/2Q p + 2kpX = irk, where k is an integer. 
Energy of zero modes in coordinate representation reads 

E = ^ E / dx J dx'v F K- 2 (x -x',n- n')d x $W(x, n)d x ,^\x', n') 

= \Y, j dx j dx'v F K- 2 {x -x',n- ri)p<®(x, n)pJV, n') (25) 

with K v is determined by Eq. (5). 

Analytical solution of Eqs. (23-24) is not a simple matter. Therefore, we consider a 
simplified model in which we neglect variations of the interaction when each coordinate varies 
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inside a segment between the impurities. That is, we approximate interaction parameter 
K~ 2 (x — x', n — n') by its value, K~ 2 (\ — j), spatially averaged with respect to coordinates 
x, n and x 1 ', n' inside given segments labelled as i and j. This model does not violate effects 
of long-range interaction which determine conductivity at low temperatures calculated in 
the next section. It is not difficult to find solution of Eqs. (23-24) for such model interaction. 
The solution for the zero-mode part of the phase operators Q u has a form 

e% ) = e v + v F '£ l K?(i-i)(8</> J 8 vp -icAN vi )t. (26) 
j 

This expression for the phase operator enables us to calculate the energy levels of zero- 
mode states given by eigenfrequencies of Green's function. This is not difficult because 
contributions from zero modes are separated from excitation in the Green's function similar 
to the case of the short range interaction considered in the Appendix. So for the energy 
levels in segment i we find 

£i = -^E K 2 (i ~ j) W„ " ttAJ^j) ± ^ E iC 2 (0) (27) 

where signs ± are, again, related to electron and hole states, respectively. 
For i = j the averaged interaction parameter can be estimated as 

A V»4( 1 + S>5)' < 28 > 

where we used dimensional units with % for clarity. 

The second term in Eq. (28) for typical values of Fermi velocity in quasi- ID conductors is 
quite large. For v F « 2 x 10 7 cm/s, which is typical value for transition metal trihalcogenides, 
its value is about 80. This corresponds to the case of strong interaction and leads to quite 
large values of coupling parameters. So we can estimate the typical energy of zero-mode 
levels as 

e 2 

Since interaction factor K p has the long-range contribution, the energy of "single-electron 
local levels" in Eq. (27) are shifted due to interaction with charges in other segments. 
Further, because of the slow decay of interaction factor K~ 2 with distance described by the 
Coulomb law, the summation in Eqs. (26,27) may diverge. In particular, (O 2 ,;) and (ef ) 
diverge for random distribution of "number of particles", y/2(5(pj/n — AiV^j) = £j — n^, at 
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segments j. The problem can be resolved by correlated distribution of localized charges in 
different segments resulting in the Coulomb gap at the Fermi energy similar to the case of 
localized shallow impurity levels in semiconductors 30 ' 31 . Indeed, Eq. (27) has some similarity 
to the energy of a local impurity level in the potential induced by all other impurities. 
Furthermore, energy of zero-mode local levels has a form somewhat resembling the expression 
for the electrostatic energy of localized impurities: 

E = T~ E(W"p " TAJ\ki)iC 2 (i - j) W,„ - TrAAy . (29) 

So we adopt to our case arguments by Efros and Shklovskii 30 for the Coulomb gap origi- 
nally used for a system of shallow impurities in semiconductors, though there are considerable 
differences between equations used in Ref. 30 and Eqs. (26-29). First, the latter contain 
operators of number of extra electrons, AN U {. But at low temperatures only a single eigen- 
value of these operators plays a role, because the states with other egenvalues n ai correspond 
to energies much larger than temperature. So we can consider the operators as c-numbers. 
Another difference is the presence of the terms related to the spin channel. Below we will 
see that this does not change the result. 

So following Ref. 30 we consider process of transfer of an electron with spin s from 
segment i in the ground state to segment j. An increase of the energy of the system induced 
by such transfer calculated from Eq. (29), with Eq. (27) taken into account, is equal to 

AE = £j - e, - — + nVFn ^ s ± 5n ^ _ KV F n ai (s - 5n ai ) > Q ^ 

where Sn a i is the difference of n a \ values in the state % after and before the transfer of the 
electron, and ry is the distance between segments i and j. Since for the process of transfer 
considered here, variation of n a \ in the initial state due to removal of an electron with spin 
s equals bn a \ = —s, while addition of the electron to the final state is described by 5n a - } = s, 
two last terms in (30) vanish, and we obtain the relation similar to that derived for the case 
of impurity levels in semiconductors. One can see that small energy differences between the 
zero-mode levels can be found only for segments situated far from each other, 

r u > (31) 

As it was shown in Ref. 30, this results in a minimum in the density of localized states, 
g(e) oc e 2 near Fermi energy, i. e., to the Coulomb gap. In particular, for electron transitions 
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with energies AE -C vpjl important at low temperatures, the distance between the segments 
participating in electron transitions must be much larger, than the typical segment length. 
For such remote segments our simplified model, in which we neglected the spatial variation 
of the interaction within distance of the order of segment length, is asymptotically exact. 

C. Effect of long-range interaction on excitation spectrum 

Solution of equation of motion (9) with long-range interaction (5) satisfying boundary 
conditions (13) is similar to the problem of interacting segments of elastic strings of random 
length. Analytical solution of this problem is difficult. So, again, we consider a simplified 
model in which the impurities are nearly arranged in planes perpendicular to the conducting 
chains, so that positions of the impurities differ by the value much smaller, than the mean 
distance between the impurities. Solution of Eqs. (8,9) consists of two parts. The first part 
corresponding to particular solution with boundary conditions (13) describes zero modes, 
their spectrum having the form of the randomly distributed localized levels discussed above. 
Levels are situated randomly because the values of 2k F Xi which determine positions of zero- 
mode levels are random. 

The part describing the excitations corresponds to general solution with zero boundary 
conditions. The latter can be found by means of Fourier transformation with respect to 
chain numbers. We take into account that according to Eq.(4) only excitations with long 
wavelengths along the chains correspond to low lying excitations (otherwise frequencies are 
of the order of plasma frequency). So we consider the limit q±_ ^> q\\ and neglect q\\ in the 
denominator of the interaction parameter K p (5). Then we easily find expressions similar 
to the last terms in Eqs. (15-16) describing the excitations, but with interaction factors and 
eigenfrequencies depending on q±, 



Thus the spectrum of the modes in this case consists of the bands, the minimum excitation 
energy being fixed by the maximum value of q± determined by the transverse reciprocal 
lattice vector ~ ir/d. Its characteristic value can be estimated as 




(32) 




(33) 
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where n = u p i/vf is Thomas- Fermi screening length in the metallic state. 

Contribution of the excitations to the Green's functions can be calculated similarly to the 
case of the contact interaction presented in the Appendix. Factors in the square brackets 
in functions B in Eq. (A4) originate from summations over eigenfrequencies in products of 
phase fields and 0„. If interaction parameter K p depends on these factors acquire a 
more complicated form, e. g., the first factor in square brackets in B must be substituted 
for 

exp (-5(^ + A - + 2 ) ln ( 1 - e ^)> qi ' (34) 

where z = n(x a — v v t a )/l. Similar substitutions occur in other factors in square brackets in 
Eq. (A4) 

In the next section we will need Green's functions for calculation of conductivity at large 
voltages, so we need to calculate Green's functions at energies e 3> uJ . Then in time and 
coordinate representation we need Green's functions for v u t a ,x a <C I. In the same time we 
consider v u t a ,x a ^> a since the energy is small compared to Fermi energy. In this limit 
expression (34) reduces to factor oc z~^ s+1 ^ with 

* = i + K * - 2 > qi (35) 

Since density of states N(e) is determined by Fourier transformation of Green's function in 
t a at x a = 0, and distribution function, n(p), of electrons on momentum along the chains 
can be found by means of Fourier transformation with respect to x a at t a = 0, this means 
that at energies e 3> lJq N(e) oc e s and n(p) oc (p — k F ) 5 , similarly to the case of Luttinger 
liquid without impurities (cf. Refs 1,3 and 11). Furthermore, if we neglect weak logarithmic 
dependence of v p in Eq. (34) on qj_, the we find that Greens function in the considered limit 
behave like that in the unbounded Luttinger liquid at any relation between v u t a and x a . In 
other words, for such relatively large energies effect of impurities on Green's function and, 
hence, on energy spectrum can be neglected. 

So we conclude that in the electronic system with Coulomb interaction the electronic 
spectrum consists of discrete levels due to zero modes, and of quasi-continuum spectrum of 
excitations of the Luttinger liquid at energies more than ZJ apart from Fermi energy. 
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D. Effect of inter-chain hopping 



Now we demonstrate that inter-chain hopping of electrons do not destroy the description 
of quasi-lD conductor with impurities at low temperatures based on the picture of collection 
of bounded Luttinger liquid. We take into account the hopping in a standard way, consid- 
ering the perturbative series for Green's functions with respect to the hopping integral t±. 
It is not difficult to calculate corrections to the Green's functions due to hopping for small 
energy e <C cJ , because at such energies Fourier-transformed Green's functions presented 
in the Appendix acquire a simple form. For example, at such energies the zero-temperature 
casual Green's function reduces to a form 



a 



[(a 2 + 4 sin 2 q xi) (a 2 + 4 sin 2 q x 2 )] " ( *<> Kp) e ikFXa+iqiXa 



G— = [V ' " ™" LJK " ' — , (36) 

2na{E — S\ + iO ■ signe) 

where e\ and q\ are energy and wave vector of the zero-level state at segment number i (c/. 
Eq. (21)). 

So calculating in a standard way corrections to the mass operator in the second order in 
t± we find that hopping results in a shift of the zero-mode level by 

&o~£^(fr (3T) 

J J 

where e- s are zero-level energies at neighboring segments. 

Note that the probability to find similar values of zero-state energies at neighboring 
segments is negligible, so the energy difference in the denominator of (37) is of the order 
of typical energy of the zero- level, uj q . So the shift of the energy due to hopping can be 
estimated as 

feo~^(T) • (38) 



LOq V I 

The higher order terms give smaller corrections proportional to higher order powers of 
factor 



f i \ 2 fa^ 



(39) 



,w J \l 

So corrections are small provided this factor is small. 

Note that our calculations are based on discrete nature of the electron spectrum near 
the Fermi energy, because we used an assumption that both characteristic energy and the 
hopping matrix element, t±, are much smaller than typical energy of zero- mode levels uJq. 
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The characteristic energy is determined by temperature, T. Thus inter-chain hopping does 
not destroy the Luttinger liquid picture in the limit T <^uJ , which does not exist in pure 
infinite Luttinger liquid. In the opposite limit of larger temperatures there are many levels in 
the energy interval of the order of temperature, and, hence, the continuum limit is applicable, 
our approach is not valid. So at temperatures T u7 , the discreteness of the excitation 
spectrum cannot be neglected, hence, according to Refs. 4,6,7,8,9, inter-chain hopping is 
expected to give significant contributions and to destroy the Luttinger liquid. 

IV. ELECTRON TRANSPORT 

In this section we calculate current in the Luttinger-liquid state of quasi-lD conductors 
by means of Keldysh diagrammatic approach 29 . 

A. Conductivity at low voltages 

Consider, first, Ohmic conductivity due to hopping between the segments confined by 
impurities. Taking into account the similarity with the case of impurities in semiconductors, 
discussed above, it is natural to expect that the mechanism of conductivity is phonon assisted 
variable-range hopping. 

There are two ways for electron transfer: via inter-chain transitions, and via transitions 
along conducting chains through potential barriers induced by impurities. In case of hopping 
conductivity there is no principal difference between both channels. So we concentrate on the 
case when probability of transitions between conducting chains dominates over probability 
of tunneling through impurities. Current flowing from segment i to its neighbors can be 
calculated as 

/ oc t ± Y^[F iti+k (l, 1) - F i+k)i (l, 1)], (40) 
k 

where -Fii + k = + Gj^ k , and summation is performed over neighboring segments. 

Green's functions Gff +k can be calculated by means of diagrams presented in Fig. 1. 

Having in mind that conductance is determined by the states with energies in vicinity of 
the Fermi energy we can use in the diagrams electron Green's functions in the form given by 
Eq. (36). Performing calculations we find that the total current through segment i consists 
of contributions describing currents flowing via segments j. Expression for such current from 
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r+t 



a 



FIG. 1: A diagram for calculation of the current. The dashed line symbolizes the phonon propaga- 
tor, Greek letters denote Keldysh time indices "+" and "— " , and a means the time index opposite 
to a. Crosses denote hopping matrix elements t±. Dots between the crosses stand for a way 
connecting segments i and j corresponding to sequential transitions via the neighboring segments. 
The lowest order contributions to current in t± are given by diagrams with the least number of 
transitions. 

segment i to segment j induced by voltage difference has a form that has very transparent 
physical interpretation 



Here and n^n- 3 are electron energies and distribution functions for zero-mode levels at 
the respective segments, and N p h(uj p h) is the distribution function of phonons. It is clear 
that different items in Eq. (41) describe processes of electron hopping between segments 
accompanied with phonon emission or absorbtion. The first factor under the integral origi- 
nates from virtual transitions between the segments along the way between segments i and 
j, and power index m describes the number of virtual transitions 



Contribution of each virtual transition to the current is reflected by factor similar to that 
presented in Eq. (39). 

Energies in arguments of the distribution functions in Eq. (41) are shifted by values 
of electro-chemical potential at given segment. Linearizing expression (41) with respect to 
potentials and performing the integration over phonon frequencies we find expression for 
the current between segments i and j similar to that describing current between localized 




+[ni(l - 7ij)(l + N ph ) - nj(l - n{)N ph ]5(ei - ej - uj ph + Viy)} . 



(41) 




(42) 
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impurity states in semiconductors 30 ' 31 



7 .. = Hi 



Rij oc exp 



2m In 



^ (L 

t±\a t 



+ 



|gj-gjl 
T 



(43) 



This expression leads to different results for the case of short-range interaction of the 
Tomonaga-Luttinger model, and for the case of long-range Coulomb interaction of elec- 
trons at different segments. In the former case, following standard arguments of theory of 
variable-range hopping 30 ' 31 we arrive at Mott's law describing three-dimensional variable- 
range hopping conductivity 



a(T) oc exp 



Tm 



1/4- 



M 



iu In 3 



0^0 / l_ 

t±\a i 



(44) 



The result is different in case of long-range interaction. Then there is the Coulomb gap 
in the density of states at the Fermi energy, and hopping of electrons occurs effectively only 
between the segments situated far from each other, at distances given by condition (31)). 
Since distances between impurities, 7, is much larger than the interchain distance, d, a given 
distance r\\ = 



|ni — nj| 2 Gp + (xi — Xj) 2 can be found for much smaller values of m, (42), 
if the segments are situated at the same or at adjacent chains, than for segments situated 
at remote chains. As current (43) between the segments exponentially decreases with m 
increasing, the strong dependence of distance between the segments on direction results in 
effectively one-dimensional hopping. Using, again, the arguments of theory of variable-range 
hopping we find that conductivity obeys Efros-Shklovskii law 



<t(T) cx exp 



ES 



T 



1/2 



T, 



ES 



ujq In 



UJQ 



I 



(45) 



Nearly ID character of hopping results in enhancement of anisotropy of conductivity. The 
anisotropy can be estimated as the ratio of the characteristic hopping lengths in directions 
parallel and perpendicular to the chains. It was argued above that hopping in the perpen- 
dicular direction occurs effectively over distances of the order of inter-chain distance d. The 
hopping length along the chains can be found as that corresponding to the minimum value 
of the exponent in Eq. (43). So for anisotropy of the conductivity we find 



<t_l d^JT ES T 



(46) 
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B. Non-linear conductivity 



Now we discuss electron transfer at high voltages and make rough estimates of non-linear 
current- volt age dependences. 

Let the voltage between m segments is equal to mV. If the voltage drop over the hopping 
length becomes larger then typical energy of the local levels, ~ uJo, then electron transitions 
between the quasi-continuum spectrum of excitations in segments become possible without 
phonon absorption. This happens at electric field E > ^/TT ES /el. Current between such 
segments can be calculated by means of diagrams similar to that presented in Fig. 1, but 
without phonon line. The diagrams contain virtual transitions via intermediate segments 
that contribute to the expression for current small factors similar to the first factor in Eq. 
(41). Further, at large voltages current is roughly proportional to factor (mV) 25+1 because of 
the power-law energy dependence of the density of states in Luttinger liquid (see discussion 
below Eq. (35)). 



Consider, first, the case when the voltage is not too large so that electron transitions between 
the regions of quasi-continuum spectrum of excitations are possible at m > 1 only. Then 
index m is related to the average electric field along the chains by relation meV ~ meEl > 
ZJq. So the most rapid, exponential, dependence of current on electric field originates from 
increase of m in the first factor in Eq. (47), in other word, from switching on new channels for 
non-linear current flow between more and more closely situated segments. The second factor 
in (47) gives slowly varying function of field E of the order of u7 25+1 , its field dependence 
can be neglected in comparison to the exponential growth due to the first factor. So for this 
regime we find 



Note that condition m > 1 means E < E \\. 

For current in direction perpendicular to conducting chains index m is related to the av- 
erage electric field by the relation meEd ^ uJq. Then for this direction the characteristic field 
can be estimated as E ± ~ T^s/ed. This is much larger value then E \\ which demonstrates 
strong anisotropy of the non-linear conductivity. 

At larger voltages, when E > E , the voltage drop is large enough to induce electron 
transitions between the neighboring segments. In this case the I-V curves can be estimated 




(47) 




(48) 
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as 

/ oc V 25+1 . (49) 

The power index here is different if tunneling along the chains via impurities is more effective 
than inter-chain tunneling. In the latter case the density of states at the ends of a segment 
is described by different power index, similarly to bounded Luttinger liquids 26 ' 27 ' 28 , so index 
28 in Eq. (49) must be substituted for (-^ l) q± - 



V. DISCUSSION 



Now we discuss relation of our theoretical results to real materials and to experimental 
data. 

Our calculations are based on the Luttinger model which does not take into account 
backscattering terms, i. e., 2k F Fourier component of the inter-electronic interaction. These 
terms are known to lead to two qualitatively different consequences (for a review see Ref. 
1). In case of attractive spin-independent interaction it results in the spin gap, while the 
backscattering terms are irrelevant if the interaction is repulsive. So, strictly speaking, our 
results can be applied for the case of Coulomb repulsion of the electrons at all wave vectors. 
On the other hand, the effective inter-electronic interaction can be considered as containing 
both repulsive Coulomb interaction and attraction due to electron-phonon coupling. It 
can happen that the Coulomb repulsion dominates at small values of wave vectors, while 
attraction induced by electron-phonon coupling dominates at wave vectors close to 2k F . 
Though such picture ignores retardation effects in the electron-phonon coupling it is quite 
fruitful in description of quasi-lD CDW-conductors 15 , electron-phonon interaction being 
described by means of attractive electron-electron coupling constant at 2k F wave vector. 
Note that the mean-field expression for the CDW gap 15 is similar to the expression for the 
spin gap obtained within the renormalization group approach with the cut-off length a in 
equation for spin gap estimated as a ~ l/k F . Therefore, we can conclude that a possible 
value of spin gap in such materials as NbSe3 and TaS3 can be of the order of few hundreds 
Kelvins. Since the spin gap implies long-range order in Q a field and expected values of 
spin gap are quite large, the fluctuations of this field at low temperatures must be small. 
The energy in the spin channel is minimized by \f2^ a = 2nn, which satisfies boundary 
conditions (13) for even values of integer m. So we can expect that in the presence of large 
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spin gap the electron system at low temperatures can be considered as spinless electrons in 
the Luttinger state, the spin degrees of freedom being frozen. In this case the results of the 
preceding sections are still valid, and only minor rather simplifications than modifications 
in the course of derivation of results are needed. 

Thus we conclude that our results can be applied to quasi-lD conductors NbSe 3 and 
TaS3 where transition from metallic to non-metallic behavior was observed. However, de- 
tailed quantitative comparison with experimental data is difficult. The first reason for this is 
some contradictory data on temperature dependence of Ohmic conductivity at low temper- 
atures. In Ref. 18 the temperature dependence corresponding to the variable-range hopping 
described by the Efros-Shklovskii law was observed, while in Refs. 16 and 17 a power-law 
dependence of conductivity on temperature was observed down to liquid helium temper- 
ature. Another difficulty for making detailed comparison is that the impurity density in 
the samples is not known. However, our results agree with general tendencies observed in 
these materials. Namely, more impure samples demonstrate transition from metallic to non- 
metallic behavior at higher temperatures, characteristic temperature in Eq. (48) increasing 
from few tens Kelvins in more perfect samples to few hundreds Kelvins in the most dirty 
samples. In a sample with characteristic temperature Ees ~ 80 K temperature dependence 
of Ohmic conductivity transforms from metallic to non-metallic behavior in the temperature 
range 100 - 200 K. Then we conclude that ZJ ~ 100 K, and since ZJ ~ e 2 /I, according to our 
approach we can estimate the mean segment length as 7 ~ 10~ 5 — 10~ 4 cm. According to Eq. 
(48) this length corresponds to electric field of the order of 10 2 — 10 3 V/cm for transition 
from Ohm's law to power law dependence. This is in order of magnitude agreement with 
the experimental data of Ref. 18. But to make convincing conclusions more theoretical and 
experimental work is needed. 

It is important to note that low temperature behavior of both NbSe 3 and TaS 3 in Refs. 
16,17,18,19 is very similar, in spite of the fact that relatively pure samples of these materials 
behave quite differently: NbSe 3 remains metallic in the CDW state, while TaS 3 below the 
Peierls transition becomes an insulator with quasiparticle density obeying the Arrhenius 
law. Whereas the main properties of quasi-lD conductors at higher temperatures are well 
understood 14 ' 15 , they demonstrate many intriguing properties at low temperatures which are 
still not explained convincingly. In particular, in TaS 3 at temperatures below 20 K a behavior 
typical for hopping conductivity was observed 32 ' 33 instead of Arrhenius law, and anomalous 
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behavior of the dielectric function interpreted as a new glassy phase 34 was detected at 
frequencies 1 4- 10 7 Hz. One can speculate that this unusual behavior can be related to the 
formation of the Luttinger phase at very low temperatures because characteristic energy 
ZJ ~ e 2 /7 is small due to small impurity density N imp = 1/7. 
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APPENDIX A: SINGLE-ELECTRON GREEN'S FUNCTION 

Here we present Green's functions introduced by Keldysh, G + ~ and G~ + . 

<G+7(U) = (iMl)<.(2)>, -<G- + (1,2) = « S (2)V>„ )S (1)>, (Al) 

with notations 1 = (ri,xi,ii), 2 = (r 2 ,X2,t2)- 

These functions are proportional to distribution function of holes and electrons, respec- 
tively. Other Green's functions can be easily derived by applying time ordering to expressions 
for G + ~ and G~ + . In particular, the casual Green's function can be found as G = G + ~ 
at ti > t 2 and G = G~ + at ti < t 2 . 

Performing the calculations in the standard way 1 ' 3 we find 

2) = -<S(1, 2) ^exp ^-* AN ^-v^ ^ exp ^(x a -v Nv t a )>j ^ 

The function related to distribution function of electrons has a form 
n-+n o\ ti(oi\I (i(5(l)i-nAN ui )(x a -v Nl/ t a )\\ ( in(x a - v Nv t a )\ 

(1 ' 2) = l) \ exp { 75 J / exp { J (A3) 

Where -8(1,2) originates from bosonic excitations of the Luttinger liquid, and the last factors 
are related to zero modes. For r\ = +1 and r 2 = — 1 this function has a form 

5(1,2) = - n« M+x) x 



V 
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I [a 2 + 4 sin 2 q Q xi] [a 2 + 4 sin 2 q x 2 ] } 



\ — e i(u v t a -qox a )-a 



I — e H.0Jut a +qox a )-a 



X 



X 



(A4) 



where x a = xi - x 2 , x s = x\ + x 2 , t a = h - t 2 , qo = ir/h, " = 77 • 

Green's functions with other values of indices r describing branches of left (right) moving 
electrons can be found from relations 

G(n = +1, xi, n = -1, x 2 ) = G(ri = +1, xi; ri = +1, -x 2 ), 
G(n = -1, xi; ri = +1, x 2 ) = G(n = +1, -x x \ r\ = +1, x 2 ), 

G(n = -1, xi; n = -1, x 2 ) = G(n = +1, -xi; n = +1, -x 2 ). 
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